function [idx,D]=findknn(X,Y,k);